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Abstract 



In this article, Temperley's bijection between spanning trees of the square grid 
on the one hand, and perfect matchings (also known as dimer coverings) of the square 
grid on the other, is extended to the setting of general planar directed (and undirected) 
graphs, where edges carry nonnegative weights that induce a weighting on the set of 
spanning trees. We show that the weighted, directed spanning trees (often called 
arborescences) of any planar graph G can be put into a one-to-one weight-preserving 
correspondence with the perfect matchings of a related planar graph 7i. 

One special case of this result is a bijection between perfect matchings of the hexag- 
onal honeycomb lattice and directed spanning trees of a triangular lattice. Another 
special case gives a correspondence between perfect matchings of the "square-octagon" 
lattice and directed weighted spanning trees on a directed weighted version of the 
cartesian lattice. 

In conjunction with results of Kenyon (1997b), our main theorem allows us to 
compute the measures of all cylinder events for random spanning trees on any (directed, 
weighted) planar graph. Conversely, in cases where the perfect matching model arises 
from a tree model, Wilson's algorithm allows us to quickly generate random samples 
of perfect matchings. 

1. Introduction 

Temperley (1972) observed that asymptotically the mxn rectangular grid has about as many 
spanning trees as the 2m x 2n rectangular grid has perfect matchings (dimer coverings). 
Soon afterwards he found a bijection between spanning trees of the mxn grid and perfect 
matchings in the (2m + 1) x (2n -|- 1) rectangular grid with a corner removed (Temperley, 
1974). The second author of the present article and, independently. Burton and Pemantle 
(1993) generalized this bijection to map spanning trees of general (undirected unweighted) 
plane graphs to perfect matchings of a related graph. Here we extend this bijection to the 
directed weighted case. 

This generalized bijection can be viewed as a way of "reducing" planar spanning tree 
systems to planar dimer systems (though not vice versa in general): for any graph whose 
spanning trees we are interested in, there is a related graph whose dimer coverings are in a 
natural one-to-one weight-preserving correspondence with the spanning trees of the original 
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graph. Thus properties of spanning trees on any planar graph can be studied by considering 
the related dimer system. However, only certain dimer systems are related to spanning 
tree systems in the aforementioned way. Two important examples are perfect matchings of 
finite subgraphs of the hexagonal honeycomb lattice (combinatorially equivalent to "lozenge" 
tilings of finite regions; see e.g. Kuperberg (1994)) and perfect matchings of finite subgraphs 
of the "square-octagon" lattice. Both of these dimer models are in bijection with weighted, 
directed spanning trees on associated graphs. 

There are a number of important applications of our bijection. Some questions about 
spanning tree models do not seem to be amenable to direct analysis, but can be approached 
if one first translates the problem into one involving the associated dimer model and then 
makes use of tools available in that context. Conversely, some problems involving dimers are 
most easily handled if one converts them into problems involving spanning trees (though this 
can be done only for a limited class of dimer models). We now describe these applications 
in greater detail. 

One example of a spanning tree property that is easy to study after reducing the problem 
to that of dimers is the computation of certain probabilities, such as the probability that 
a directed edge ei is in the tree and the directed dual edge is in the dual tree. (For a 
definition of dual tree, see § ||.) The presence or absence of the dual edge in the dual 
tree is not a local event with respect to the (primal) tree model; that is, the event is not 
determined by the presence or absence of a fixed set of edges in the primal tree. (The fact 
that 62 is an oriented edge is crucial here.) On the other hand, the event in question is 
a local event in the associated matching process, since the matching directly incorporates 
both primal and dual directed trees. The probabilities of local events in either the tree or 
matching model are easy to compute (Burton and Pemantle (1993), Kenyon (1997b)), but 
events of the above type are harder if not impossible to compute from the point of view of 
the tree only (Burton and Pemantle, 1993). 

Another spanning tree property that can be studied via dimers is the number of times 
that the path connecting two points in a spanning tree winds around the two points. In § ^ 
we relate these winding numbers to height functions in the dimer model; the first author 
has shown in Kenyon (1997b) how to compute properties of these height functions (and 
consequently the corresponding winding numbers) such as the variance. 

Dimer systems can also be studied via trees, if the given dimer system has a spanning tree 
model associated with it. For instance, one can sometimes enumerate the dimer coverings of a 
graph by counting the number of spanning trees in the associated tree model. In § ^ we show 
a variety of such graphs, together with exact formulas for the number of dimer coverings, 
where the easiest (or only) way we know to obtain these formulas is to count spanning trees. 
In the dimer model on a bounded region, the boundary can have an important (long-range) 
effect on the number of configurations (Cohn et al., 1998). In this case the regions which arise 
from the associated spanning tree process give the most "natural" boundary conditions for 
the dimer model, in the sense that the boundary has the least long-range influence (Kenyon, 
1997a). 

Another case where a spanning tree model is useful for studying the associated dimer 
model is in the generation of random samples. Wilson's algorithm (Propp and Wilson, 1998) 
can be used to generate random spanning trees quickly — the expected running time is given 
by the sum of two mean hitting times. For the lattice of octagons and squares, the expected 
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running time is linear in the number of vertices. For the usual lattice of squares, when a 
rectangular region has moderate aspect ratio, the running time is nearly linear, but with a 
logarithmic correction factor. 

Finally, Burton and Pemantle (1993) prove that the uniform measure on spanning trees 
of the n X n square grid converges as n ^ oo to the unique translation-invariant measure 
of maximal entropy on the set of spanning forests with no finite component. Consequently 
the associated dimer model (on Z^) has a unique translation-invariant measure of maximal 
entropy. We do not know how to prove this directly from the dimer model itself, or in any 
other dimer model except those arising from our construction via a bijection with undirected 
(but possibly weighted) spanning trees. 

We remark that other combinatorial systems that can also be reduced to dimer systems in 
a similarly simple way include the Ising model on planar graphs (Fisher, 1966) and systems 
of non-intersecting lattice paths (Lindstrom, 1973; Gessel and Viennot, 1989). 

In § we prove the generalized version of Temperley's bijection. In the two succeeding 
sections (§§ Q-]^) we illustrate the bijection with two examples: In § ^ we exhibit a bijec- 
tion between directed spanning trees on the triangular lattice and perfect matchings of the 
hexagonal honeycomb lattice. Our bijection cannot be applied directly to matchings on the 
square-octagon lattice, but in § ^ we show how to locally transform this lattice so that the 
bijection can be applied. This transformation enables the rapid generation of random dimer 
configurations of the square-octagon lattice. Then in § |^ we show how the winding number 
of arcs in a spanning tree can be related to the height function on the corresponding perfect 
matching. In § ^ we use our generalized bijection to compute the exact number of perfect 
matchings of some "locally symmetric" finite planar graphs, that is, graphs that arise as 
finite induced subgraphs of highly symmetric infinite planar graphs. Lastly, in § |^ we give 
some open problems. 

2. Generalized Temper ley Bijection 

Let G be a finite connected directed graph embedded in the plane, with multiple edges and 
self-loops allowed. In general the edges of G will be weighted, that is, each directed edge 
from vertex u to vertex v has a nonnegative weight assigned to it, which need not be the 
same as the weight of other directed edges from m to or from v to u. Undirected graphs 
can be fit into our framework by thinking of each undirected edge as two directed edges, one 
in each direction, embedded in the plane so as to coincide. (We will discuss issues related 
to choice of embedding at the end of this section.) Unweighted graphs can be fit into our 
framework by assigning each edge weight 1. 

By a directed spanning tree (or arborescence) T of G we mean a connected, con- 
tractible union of (directed) edges such that each vertex of G except one has exactly one 
outgoing edge in T. Note that the exceptional vertex necessarily has no outgoing edges in 
T; it is called the root of T. We define the weight of such a tree T to be the product of 
the weights of its edges. 

We will make a new weighted graph 7i(G) based on G, as shown in the top half of 
Figure G is shown in the top left-most panel. G^, the dual graph of G (second panel), has 
vertices, edges, and faces of G^ corresponding to faces, edges, and vertices of G, respectively 
(including a vertex, here marked /*, that corresponds to the unbounded, external face of 
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Figure 1: Illustration of the generalized Temper ley bijection. 
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G, and is represented in "extended form", spread-out region rather than a small 

dot). We can embed G and G-^ simultaneously in the plane, such that an edge e of G crosses 
the corresponding dual edge e"*" of G^ exactly once and crosses no other edge of G"*". If we 
introduce a new vertex at each such crossing, we get the graph shown in the third panel. 
This is the graph 7i(G). Pictorially, we may derive 'H{G) from G by adding a new node on 
each edge e and a new node on each face / and joining them by a new edge if e is part of 
the boundary of /. To avoid confusion, we will say that 7i(G) has nodes and links whereas 
G has vertices and edges. 

Here is an alternative, direct definition of 1-L{G) that does not go by way of the dual 
graph. Put V = the set of vertices of G, E = the set of edges, F = the set of faces 
(including the unbounded face). Define H{G) as the weighted undirected graph with a node 
V corresponding to each vertex v of G, a node e corresponding to each edge e of G, and 
a node / corresponding to each face / of G, with a link joining two nodes in T-C{G) if the 
corresponding structures in G are either an edge and one of its endpoints or an edge and 
one of the faces it bounds. 

The weight of a link between a vertex-node v and an edge-node e (where v is an endpoint 
of e in G) is the weight of edge e in G directed away from v. The weight of a link between 
a face- node / and an edge-node e (where e bounds face / in G) is always 1. 

A perfect matching of a graph if is a collection of edges M such that each vertex is 
a vertex of exactly one edge of M. The weight of a perfect matching is the product of the 
weights of its edges (1 by default in the unweighted case). 

In the case of both trees and matchings, the weighting gives rise to a probability distri- 
bution on the objects in question, in which the probability of any particular object (tree or 
matching) is proportional to its weight. 

Let V* be a vertex of G and /* a face of G, and let H = Tiiy* , f*) be the induced 
subgraph of TiiG) obtained by deleting the nodes t>*, /* (along with all incident edges in 
T-l{G)), as shown in the fourth panel of the top half of Figure |l|). Since by Euler's formula 
{\V\ — 1) + {\F\ — 1) = H{v*, /*) is a balanced bipartite graph, so it may have perfect 
matchings. (For a nice tree-based proof of Euler's formula, see (Aigner and Ziegler, 1998, 
page 57).) 

Theorem 1 If v* is incident with f* , then there is a weight-preserving bijection between 
spanning trees of G rooted at v* and perfect matchings of T-C{v* , f*) . If v* is not incident 
with f*, there remains a weight-preserving injection from the spanning trees of G rooted at 
V* to the perfect matchings of Ti^v* , f*) . 

This theorem, along with its proof, is a generalization of a result of Temperley (1974) 
which is discussed in problem 4.30 of (Lovasz, 1979, pages 34, 104, 243-244). The unweighted 
undirected generalization was independently discovered by Burton and Pemantle (1993), who 
applied it to infinite graphs, and also by F. Y. Wu, who included it in lecture notes for a 
course. 

Note that in the special case when we take all weights of G to be 1, the first part of the 
theorem implies that the number of perfect matchings of 7i(f *, /*) is independent of v* and 
/*, provided that v* and /* are incident with one another. 

Henceforth, we refer to perfect matchings as simply "matchings," and directed spanning 
trees of G rooted at v* as simply "spanning trees" or occasionally just "trees." 
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Proof of theorem: It will be enough to exhibit a weight-preserving injective mapping 
from the set of spanning trees of G into the set of matchings of T-C{v*, f*), and to show that 
when V* is incident with /*, every matching of Ti.{v*, /*) arises from a spanning tree of G. 

Given a spanning tree T of G rooted at V*, the set of edges of G"*" that do not cross edges 
of T form a spanning tree of G^, called the dual tree and here denoted by T"*". Orient the 
edges of T-*- so that they point towards /*. Then a matching M of Ti.{v*, /*) can be obtained 
as shown in the bottom half of Figure |l|. Specifically, for each v E V, pair v with the unique 
e such that v is an endpoint of e and e is pointing away from v in the orientation of T, and 
for each f E F, pair / with the unique e such that e bounds / and e-*- is pointing away from 
/ in the orientation of T-*-. The left panel shows the tree T; the second panel shows the 
dual tree T-*-; the third panel shows both trees; and the fourth panel shows the matching M, 
which has the same weight as T. 

To verify that this construction always gives a matching M of 7Y(f *,/*), it suffices to 
show that no edge-node e is paired twice. But this could only happen if we had e G T and 
g± e T"*-, contradicting the definition of a dual tree. 

From the matching M we can easily recover T as the set of edges e such that e is paired 
with a vertex-node in Ti.{v*, /*) under the matching M. Hence the mapping T \—>- M is 
injective. 

Now suppose V* is incident with /*, and let M be a matching of 7i(f*, /*). Let T be the 
set of edges e of G such that e is paired with a vertex-node by M. To complete the proof of 
the theorem, we must show that T is a spanning tree. Note that T has — 1 edges, so it 
suffices to prove that T is acyclic. 

Suppose T contained a cycle C, say of length n. C divides the plane into two (open) 
regions, one of which contains both v* and /* and the other of which contains neither. We 
claim that each part contains an odd number of nodes of 7i(G) and hence an odd number 
of nodes of the subgraph H{v*,f*) as well. For, suppose we modify G by replacing either 
of the two regions by a single face. By Euler's formula, the number of vertices, edges, and 
faces in the resulting graph must be even. Since there are an even number of these elements 
on the cycle C {n vertices and n edges) and an odd number in the modified region (1 face), 
the unmodified region must have an odd number of elements as well. 

Since the edges of C disconnect H{v*,f*) into parts lying in the two regions, M must 
match each region within itself. But this is impossible, since each region has been shown to 
contain an odd number of nodes of 7i(f *, /*). This completes the proof of the theorem. □ 

As was remarked earlier, the theorem implies that when v' is incident with /' and v" is 
incident with /", the matchings M' of 7i(t'', /') are equinumerous with the matchings M" 
of 7^(f", /"); in fact, the proof of the theorem provides a bijection between the two sets of 
matchings. This bijection can be understood without reference to spanning trees, as a process 
of "sliding edges." Specifically, one iteratively defines a chain v" = vq, cq, fi, ei, V2, ... such 
that, for all i, is the node that M' pairs with and f j+i is the vertex of Cj that is distinct 
from Vi. This chain cannot repeat any vertices, since any closed loop would encircle an odd 
number of nodes (see the preceding proof), so it must terminate by arriving at v' after some 
number of steps. That is, the chain must be of the form 



v" = Wo, eo, wi, ei, V2, Cr-i, v,. = v' 
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for some r. Once one has found such a chain, one modifies the matching M' by pairing 
with instead of Vi. One then does the same with a chain of dual-edges joining the faces 
/" and /', obtaining the desired matching M". 

We also remark that in addition to one's having a choice of which vertex-node and face- 
node to delete, one often has a choice of how to embed a graph in the plane in the first place. 
For instance, in the case where G has a single edge from m to f and a single edge from v to 
u, we allowed the two edges to be embedded so as to coincide. What if we had required the 
embedding to be proper, so that the two edges could meet only at their endpoints? Then 
one would get a slightly enlarged graph T-C{G) in which a single edge-node in the original 
7Y((j') was replaced by two edge-nodes and a face-node in between (corresponding to the 
digon bounded by the two edges). It is easy to see in this case that perfect matchings 
of the first T-C{G;v*, f*) are in bijection with perfect matchings of the second TC{G;v*,f*). 
When there are multiple directed edges in each direction, the number of possible embeddings 
increases rapidly; but our main bijection theorem guarantees that the number of matchings 
of H{G; V*, /*) is insensitive to the choice of embedding. 

Moreover, having several directed edges from f to w is in a certain sense equivalent to 
having a single edge from v to w whose weight is the sum of the weights of those directed 
edges. It is not true that the spanning trees of the former graph are in bijection with those 
of the latter graph; however, there is an obvious mapping from the former to the latter, and 
this correspondence is weight-preserving, in the sense that the weight of a spanning tree of 
the smaller graph is the sum of the weights of the spanning trees in the larger graph to which 
it corresponds. It follows that the sum of the weights of all the spanning trees is the same 
for both graphs. 

Given a graph H, it can be an amusing problem to find a directed graph G such that 
Tt{G;v*,f*) = H. We leave it to the reader to show that this cannot be done with the 
square-octagon lattice of Figure |^. (That is, there is a finite subgraph of the lattice, such 
that any subregion H of the square-octagon lattice containing this subgraph will fail to be 
of the hTmniG;v*,f*).) 

3. The Hexagonal Lattice 

In this section we illustrate the technique of Theorem |1] by giving a bijection between span- 
ning trees of a directed graph and matchings in the hexagonal (honeycomb) lattice. 

Panel (a) of Figure ^ contains the plane graph G, a directed triangular lattice. (Here 
and throughout the rest of the article, the upper-left, upper-right, lower- left, and lower-right 
panels of a four-panel figure will be denoted by (a), (b), (c), and (d), respectively.) G 
contains an "outer vertex" which is represented in extended form, in this case drawn as a 
large hexagon. In the examples throughout the rest of the article, either G or G-^ (or both) 
will have an outer vertex that is drawn in extended form. Panel (b) shows the dual of G, a 
hexagonal lattice. The edges in panel (b) have been drawn bent slightly so that the union 
of G and its dual (panel (c)) can be recognized as a subset of the hexagonal lattice. The 
dotted edges in panel (c) have weight zero, and may be omitted; they are shown only to 
highlight the connection with panel (a). The graph 'H{G) can be read off from panel (c); it 
is a hexagonal lattice with about three times as many hexagons as G^. Panel (d) shows the 
graph H{v*, /*), which is obtained from 7i(G) by removing v* (the outer vertex of G) and 
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Figure 2: Generalized Temperley bijection for the hexagonal lattice. 
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/* (the leftmost vertex at the top of G^). We shall apply this correspondence in 



4. The Square- Oct agon Lattice 

Here we illustrate a less direct application of Theorem |I], and give a bijection between perfect 
matchings of certain planar graphs and spanning trees on an associated graph. Consider 
perfect matchings on the square-octagon lattice, an excerpt of which is shown in Figure ^ 
This graph does not arise as 7i(G) for any graph G, so Theorem |l] does not apply immediately. 
Nonetheless, it is possible to generalize the bijection to apply to this lattice. To do it 



Figure 3: A portion of the square-octagon lattice. 



we need to apply two transformations to the lattice. The first transformation is called 
"urban renewal," a term coined by the second author, who learned of the method from Greg 
Kuperberg. In the second transformation, we adjust the edge weights. At that point, if 
a square-octagon region has suitable boundary conditions, the transformed graph can be 
expressed as T-C{G) for some graph G. 



4.1. Urban renewal 

Tricks such as urban renewal have been used by researchers in the statistical mechanics 
literature for decades, but since understanding it is essential for what follows, a description 
is included here of the special case of urban renewal that we will need. One views the square- 
octagon lattice as a set of cities (the squares) that communicate with one another via the 
edges that separate octagons. Now the graph of cities (with each city being thought of as 
adjacent to the four closest cities) is itself bipartite, so we may say that every city is either 
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rich or poor, with every poor city having four rich neighbors and vice versa. The process 
of urban renewal on a poor city merges each of its four vertices with its four neighboring 
vertices, and then changes the weights of the edges of the poor city from 1 to 1/2, as shown 
in Figure |[ We will show that the sum of the weights of the matchings in the "before" 
graph is twice the sum of the weights of the matchings in the "after" graph. We will do this 
by associating with each matching in the before graph one or two matchings in the after 
graph, and vice versa. More precisely, we divide the set of matchings in the before graph 
into equivalence classes of size 1 or 2, and likewise with the set of matchings of the after 
graph, and we create a bijection between these equivalence classes so that the weight of each 
class in the before graph (that is, the sum of the weights of the matchings that constitute 
that class) is twice the weight of the associated class in the after graph. 



1/2 



1/2 



1/2 



Figure 4: Urban renewal. The poor city is the inner square in the left graph, and is connected 
to the rest of the graph (not shown) via only four vertices at the corners, some of which may 
actually be absent. The city and its connections are replaced with weight 1/2 edges, shown 
as dashed lines. All other edges have weight 1. 



Matchings in the "before" graph get mapped via urban renewal to matchings in the 
"after" graph by deleting the four vertices of the poor city and its incident edges, and then 
pairing up any resulting unpaired vertices. Prior to urban renewal, every matching will 
match k of the poor city's vertices with the rest of the graph, with k equal to 0, 2, or 4; 
if = 2, then these vertices are adjacent. U k = 0, then since the city has two possible 
matchings, a pair of matchings in the "before" graph get mapped to one matching (of half 
their combined weight) in the "after" graph, li k = 2 (two of the poor city's vertices match 
to each other and two match outward), then the matching in the before graph gets mapped 
to a matching in the after graph that uses one weight-1/2 edge. The matchings with k = 4 
get mapped to a pair of matchings in the after graph, each using two weight-1/2 edges. Thus 
urban renewal on a poor city will reduce the weighted sum of matchings by a factor of 1/2. 
(If one is trying to generate random matchings rather than merely count them, then, given 
a random bit, a random matching in the before graph is readily transformed into a random 
matching in the after graph, and conversely, given a random bit, a random matching in the 
after graph is readily transformed into a random matching in the before graph.) 

The preceding discussion applies to cities in the interior of a finite subgraph of the 
infinite square-octagon grid. Along the boundaries, some of the poor cities may not have 
four neighbors, but urban renewal can still be done. One way to see this is to adjoin a pair 
of connected vertices to the graph for each missing poor city's neighbor, and connect one of 
these vertices to the poor city. This operation won't affect the number of matchings or their 
weights, and after urban renewal, the pair may be deleted again, again without affecting the 
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matchings — so if some of the poor city's vertices don't have neighbors, these vertices are 
deleted by urban renewal. 

Doing urban renewal on each of the poor cities in the square-octagon lattice will yield 
the more familiar Cartesian lattice. 

4.2. Weighted directed spanning trees 

Consider the finite square-octagon graph shown in Figure ^. It has 3 octagonal bumps on 
the left, and four on top, so by convention let's call it a region of order 3,4. (In a region 
of order L,M, there are 2LM octagons.) An octagonal column and octagonal row meet 
at a unique square; these will be the rich cities. The rich cities have been labeled by their 
coordinates to enhance clarity. The other [L + ^){M + 1) squares will be the poor cities, 
and we will do urban renewal on them as shown in Figure ^. We will compute the weighted 
number of matchings of the resulting graph, and multiply by 2(-^+^)(*^+^). 

Now for any vertex, we may re-weight all the edges incident to the vertex (multiplying 
them all by the same constant) without affecting the probability distribution on matchings: 
this has the effect of multiplying the weight of any matching by that same constant. Re- 
weight the edges as follows: for the rich city (complete or incomplete) with coordinates i,j 
(0 < z < M, < j < N), multiply the weights of the edges incident to the top left and lower 
right corners by 2^*"-', the other two corners by 2*+-'. 

Edges that are internal to the rich cities remain weighted by 2^~^^2~^~^ = 1. The long 
edges come in pairs. The lower or right edge of the pair gets its weight doubled, to become 
1, while the upper or left edges of the pair gets its weight halved to become 1/4. 

The next thing we need to do is interpret this graph as a plane graph and its dual (see 
Figure The upper left vertices of the small squares represent vertices, the lower right 
vertices represent faces, and the other two vertices represent edges. The result is the graph 
F shown in Figure |^, which has LM + 1 vertices — LM of them on a grid, and one "outer 
vertex" (not in the original graph) that all the open edges connect to. A random spanning 
tree on the vertices of this graph rooted at the outer vertex determines a dual tree on the 
faces of this graph, rooted at the upper left face, and the two together determine a matching 
of the graph in Figure |^. The weight of the matching equals the weight of the primal tree, 
since the re-weighting left every dual tree with weight one. 

4.3. Random generation in linear expected time 

Using the loop-erased random- walk spanning tree generator (Propp and Wilson, 1998), and 
the bijection derived above between spanning trees of a weighted graph and matchings of 
the square-octagon regions, we can sample random matchings in linear time. The tree 
generator builds the tree by doing a sequence of loop-erased random walks on the underlying 
graph. (From any vertex v, the probability of moving to a particular neighboring vertex w 
is proportional to the weight of the edge from v to w; this determines a random walk on the 
graph. For details on loop-erasure, see (Propp and Wilson, 1998).) It has been shown that 
the expected running time (or rather, number of random-walk steps) of the tree algorithm 
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Figure 5: Region of order 3, 4 before and after urban renewal. The poor cities on which 
urban renewal is done are labeled with "ur," the rich cities are labeled with their coordinates. 
Dashed edges have weight 1/2. 
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Figure 6: The region from Figure ^ after re-weighting. Dashed edges have weight 1/4. A 
distinguished vertex and distinguished outer face have been adjoined to give a graph H which 
nicely decomposes into a graph G and a weighted directed graph F. 



is given precisely by 



# times a that random walk started 
at V visits v before hitting the root 



For our random walk, the moves are right or down with probability 4/10 and up or left 
with probability 1/10, since in the face graph the links going to the left or up have 1/4 the 
weight of the other links. For large graphs, the random walk drifts to the right and down, 
so we consider this biased random walk on Z^. Starting at the origin, with probability 1 the 
origin is visited finitely many times. Let R be the expected number of times the random 
walk returns to its starting location, counting the "return" at time 0, before drifting off to 
infinity. The first expression below for R is not hard to check, and the remaining two can 
be found in (Beyer, 1981, p. 408): 

2 



fc=0 

»1 



(1/5)^' 



JO 



dx dy 



1 - (2/5) cos(7rx) - (2/5) cos(7ry) 
= (2/7r)ir(4/5) 

with K{k) denoting Legendre's complete elliptic integral of the first kind. The expected 
number of steps to create the random spanning tree is bounded hy R = 1.27025 times the 
number of vertices, and the remaining steps are readily done deterministically in linear time. 
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5. Height Functions and Winding Numbers 

In this section we describe the connection between the winding number of a spanning tree 
on a planar graph G and a height function on its corresponding matching graph Ti. The 
result, Theorem |^ below, answers a question posed to the first author by Itai Benjamini. 

5.1. Height function definition 

We assume that G is connected and that T-C{G) is embedded in the plane with straight 
edges. The straight-line embedding is not necessary for the definition (see below) but the 
construction is more geometric in this case. Moreover, we assume that T-C{G) is embedded 
so that one of /* or v* is the outer node (as in Figure |^), or else both /* and v* are on the 
outer facet (as in Figure ||). 

Recall that each facet of T-C{G) is a quadrilateral containing a node v, a node /, and 
two nodes 61,62. The nodes v and / are opposite each other. Let d be the diagonal of the 
quadrilateral facet directed from v to /. Let arg((i) G [0, 27r) denote the angle of the vector 
d with respect to the x-axis. 

Let M be a perfect matching of 'Hiy* , /*) and T, T-*- respectively the associated spanning 
tree and its dual. Let D be the set of the diagonals of facets of T-C{G). We will define a 
real- valued height function : D ^ M associated with matching M (refer to Figure |^ and 
Figure ^ . 

Remark In many of our examples in § |^, one or more vertex or face nodes are drawn in 
an extended format. In these cases, the "diagonals" incident to an extended node may be 
drawn from any point in the node. In many situations it is natural to draw more than one 
diagonal on a facet if one of the nodes bounding it is drawn in an extended fashion (as in 
Figure |^), and then each of the diagonals gets its own height. For instance, to recover the 
standard definition of height function for matchings of subgraphs of 7?, it is necessary to 
draw multiple diagonals (see Figure ||). It is thus more natural to view the heights as being 
defined on the diagonals of the facets rather than the facets themselves. 

We first cut the plane along the links in the perfect matching M. We need every vertex- 
node and face- node to be at the end of one cut, and since neither v* nor /* is in the matching 
M, we make one additional cut, from v* to /*. If both v* and /* border the outer facet, 
then we make the cut between them pass through 00. If one of v* or /* is the outer node, 
then we view the outer node as being at 00, so that the cut still goes to 00. It is convenient 
to make this cut split the diagonal from v* to /*, so that there are two diagonals in D from 
V* to /*, one on either side of the cut. 

We require the height function h to satisfy the following local constraint. Suppose d and 
d' are two diagonals that share a vertex node or face node x. Let 6 be the angle required 
to rotate d to d' around x. Either the counterclockwise (positive) rotation or the clockwise 
(negative) rotation will encounter the cut containing x; we take 6 to be the rotation which 
avoids the cut. The local constraint is 



h{d') = h{d) + e. 
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Figure 7: The height function for the perfect matching in Figure |T[ Since v* and /* are 
unmatched, the height drops by 27r on the facet containing v* and /*. 




Figure 8: The height function associated with a perfect matching of another graph. For 
backwards compatibihty with previous definitions of the height function associated with 
dimers on the square lattice, we have (1) measured the heights in quarter-turns rather than 
in radians, which introduces a scale factor of 7r/2, and (2) drawn multiple diagonals (each 
with its own height) on facets bounded by an extended node. 
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Lemma 2 Up to a global additive constant, there is a unique height function which satisfies 
the local constraints. The additive constant can be chosen so that for each diagonal d, h{d) = 
arg((i) mod 27c. 

Proof: If two diagonals share the same vertex node or face node, then their height 
difference is determined by the local constraints. Suppose that vertices Vi and f 2 of G are 
connected by an edge, and that / is a face bounded by this edge. The height difference 
between the diagonal from vi to / and the diagonal from V2 to / is determined, and conse- 
quently, the local constraints determine the height difference between any diagonal with vi 
as vertex node and any other diagonal with f 2 as vertex node. Since G is connected, there 
can be at most one height function (up to global additive constant). 

We next check that the local constraints do not overconstrain the height function, i.e. 
that there is a height function satisfying them. Between any two diagonals for which there is 
a local constraint, one can draw a path connecting the diagonals but which avoids the cuts. 
If the local constraints were inconsistent, then there would be a closed loop in the plane, 
which avoids the cuts, such that the local constraints on the diagonals crossed by the loop 
are inconsistent. Since the cut from v* to /* passes through 00, the interior of this loop does 
not contain v* or /*. Consider such a contradictory loop surrounding a minimal number of 
links in the matching M. The loop must cross at least one diagonal between a vertex node 
and a face node, and one of these nodes must be in the interior of the loop. Call this node 
X. Since x is not v* or /*, it is paired with an edge node y in the matching. Since the loop 
avoids cuts, y is also in the interior of the loop. 

Since x is a vertex node or a face node, the local constraints on the diagonals incident to 
X involve rotations that avoid the cut from x to y, and they are evidently consistent. Since y 
is an edge node, there are four facets of H{G) incident to y, and the four diagonals in these 
facets form a quadrilateral containing y. The total height change going clockwise around y 
is then the sum of the interior angles of this quadrilateral, excluding the angle at the node x. 
Since the sum of the interior angles of a quadrilateral is 27r, the total height change around 
y is 27r minus the angle at x. Thus the total height change going around the cut from x to 
y is constrained to be zero. Therefore the contradictory loop can be deformed to exclude x 
and y from its interior, contradicting are assumption that it surrounds a minimal number of 
links from the matching M. We conclude that there are no such contradictory loops, and 
that the height function is well-defined (up to a global additive constant). 

The second statement of the lemma follows by noting that if for some diagonal d, h{d) = 
arg((i) mod 27r, then this relation holds for the neighboring diagonals as well. □ 

Note that in the case of a matching of Z^, this definition of height function is 27r/4 times 
the standard definition due to Thurston (1990) (see Figure H). It is also essentially equivalent 
to, but more geometrical than, the definition due to Propp (1993). 

When Hi^G) is embedded in the plane but not with straight edges, one can still assign 
to each diagonal d an angle which is the argument of the vector representing the difference 
in its endpoints. The orientation of the triple d, i, d' is a topological condition and so does 
not depend on the fact that the edges of 'H(G) are straight. Thus it is possible to define the 
height for any embedding of 1-L{G). 
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5.2. Turning and heights 

Let 7 be a simple path (in topological terms, a "directed arc") in the spanning tree T from 
v' to f specifically 7 is the link path 

V' = fo, Co, Vi,ei,. . . , fr._i, Cr-l, Vr = v" . 

Note that the matching of ?i(f *, /*) that corresponds to T matches node Vi with node Cj for 
all < z < r — 1. Let be an edge- node (adjacent to v^) towards which the path can be 
continued. 

Let /i, . . . /fc be the chain of facets of T-C{G) which share a node with 7 and lie to the left of 
7, where /i contains the first link vqCo and fk contains the link e^-ifr- Then for i G [1, A; — 1] 
the facets /« and /j+i are adjacent along a single link of 7^(G), which furthermore is an 
unmatched link (this is a link which shares a node with 7 but is not in 7; therefore it is 
unmatched). 

The winding number of 7 is defined to be the total angle of the left turns minus the total 
angle of right turns, from vq to Vr (the "initial" angle of 7 is the direction of VqCq and the 
"final" angle of 7 is the direction of VrCr)- 

Theorem 3 The winding number of'j is equal to {h{fk) — c^) — {h{fi) — Ci), where Ci is the 
counterclockwise angle from the vector vqCq to the diagonal di, and Ck is the counterclockwise 
angle from the vector v^Cr to the diagonal d^- 

Here we may view the heights as being defined on the facets, since even if facet fj has 
more than one diagonal, h{fj) — Cj has the same value no matter which diagonal is used. 

Proof of theorem: Without loss of generality we may assume that the height at facet 
fi is h{fi) = arg((ii). Then h{fi) — ci is the angle that the initial direction of 7 makes with 
the X-axis. Let fj be a facet adjacent to ve. Then similarly h{fj) — Cj is equal modulo 27i to 
the angle that ViCi makes with the x-axis. 

If fj is the last facet adjacent to vi, so that fj+i is adjacent to f^+i, then the difference 
(/i(/j_i_i) — Cj+i) — {h{fj) — Cj) equals the increase in angle from VjCj to fj+iCj+i (which is 
negative at a right turn). The proof follows. □ 

6. Explicit Formulas 

Here we use the generalized Temperley bijection to count perfect matchings of certain finite 
subgraphs of the infinite square grid and infinite square-octagon grid, making use of spanning 
trees. The enumeration techniques closely follow the derivation of the exact formula for the 
number of domino tilings of the {2n + 1) x (2n+ 1) square with a corner removed; see Lovasz 
(1979) and Propp (1995). 

One motivation for some of these calculations is that they can be used to compute asymp- 
totic formulas for the number of dimer configurations on more general regions, via techniques 
developed in Kenyon (1998). For example, using the exact formula for the triangular or dia- 
mond regions in §§6.5-6.8, one should be able to extend the asymptotic formula in Kenyon 
(1998) to polygonal regions whose boundary edges have slopes in {0, 1, —1, 00}. Using the 
example in |6]^ should give rise to a similar asymptotic formula for regions of the hexagonal 
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lattice. Other formulas corroborate refinements of the entropy formula that predict how the 

the lower-order asymptotics of the number of spanning trees should reflect the geometry of 
the boundary of the graph (see Duplantier and David (1988) and Kenyon (1998)). 

To enumerate the spanning trees of a graph, we make use of the well-known Matrix Tree 
Theorem (see e.g. Biggs (1993)), as illustrated below. Given a directed graph G with n 
vertices, the negative Laplacian of G is the n x n matrix L{G), where L{G)j},w (v 7^ w) 
equals the negative of the weight of the edge from vertex v to vertex w, and L{G)y^v is the 
weighted sum of the arcs emanating from v. The determinant of the submatrix obtained by 
deleting row r and column r from L[G) gives the (weighted) number of the spanning trees 
rooted at r. 
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= 17716630 



We will make use of two ways to evaluate this determinant. We can exhibit n — 1 
orthogonal nonzero eigenvectors of the matrix obtained from L{G) by deleting row and 
column r, and multiply their eigenvalues. In the case where the graph is undirected, so that 
the matrix is symmetric, an alternate procedure is available: we can exhibit n orthogonal 
nonzero eigenvectors of L{G), multiply their eigenvalues except the zero eigenvalue, and 
divide by n. (A more accurate general description of this procedure is that the multiplicity of 
the eigenvalue zero is to be reduced by 1. When zero is a multiple eigenvalue, the procedure 
gives rise to the product 0. However, in all our examples, zero is a simple eigenvalue, 
corresponding to the fact that the graph is connected and hence possesses one or more 
spanning trees.) This second method can be viewed as a variation on the first method, in 
which an auxiliary vertex is added to the graph, the edges from every other vertex to it are 
given weight e, the trees rooted at the auxiliary vertex are counted, e is sent to 0, and the 
result is divided by n, corresponding to the fact that any of the n vertices of the original 
graph could become the unique vertex joined to the extra vertex in a spanning tree in the 
new graph. (The details are left to the reader.) 

Every graph considered in the examples below either is a finite induced subgraph of 
an infinite square grid or else is obtained from such a graph by adding a single additional 
vertex. In each case we give an explicit formulas for the eigenfunctions. We hasten to say 
that for a general planar graph an explicit diagonalization would be much more difficult. The 
tractability of our chosen examples arises from their connection with the negative Laplacian 
of the infinite square grid Z^. A function / on is an eigenfunction of the negative Laplacian 
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if it satisfies the equation 

4/ (a;, y) - f{x - 1, y) - f{x + l,y)- f{x, y-l)- f{x, y + l) = A/(x, y) 

for all integers x, y, in which case A is the associated eigenvalue. For each pair of complex 
numbers satisfying \(\ = \('\ = 1, we can construct such an eigenfunction by putting 
f{x,y) = C^C'^ and X = A — ( — ('-^ — C ~ C'~^- many cases, the restriction of such a 
function f{x,y) to points x,y lying in some finite region in is an eigenfunction of the 
matrix associated with that region by the Matrix Tree Theorem. 

As a preparatory example, we briefly mention here the trivial case of counting spanning 
trees of an undirected chain C consisting of n vertices and n — 1 edges. The space of 
complex-valued functions on the vertices of C can be identified with the ra-dimensional space 
W of odd periodic functions of period 2n + 2, i.e. functions / : Z i-^ C that satisfy f{—x) = 
f{2n+2—x) = —f{x) for all x in Z (and that consequently satisfy /(O) = f{n+l) = ■ ■ ■ = 0). 
Under this identification, the negative Laplacian of C is carried over to the negative Laplacian 
of Z, giving rise to an eigenbasis for W of the form fk{x) = sin (1 < A; < n). 

We draw each graph (see Figures ^ through |1^ so that the non-root vertices are located 
at points of the lattice, and so that every edge of the graph connects nearest neighbors in the 
lattice. With high-degree root vertices it is necessary to draw the root vertex in an extended 
fashion, covering multiple points of the lattice, to ensure that all the edges incident to the 
root are drawn between neighboring points in the lattice. 

Generally the requisite eigenvectors of the matrix L{G) (either with or without the root 
row and column deleted) will be eigenvectors of the negative Laplacian of the infinite lattice 
that also satisfy certain boundary conditions (analogous to the conditions /(O) = f{n+l) = 
in the preparatory example). Usually the infinite lattice will be Z^ although sometimes it 
will be (Z + |)^ = {(x, y) : X — |, ?/ — I G Z}. There are two types of boundary conditions. 
Suppose that v and w are nearest neighbors in the lattice, and f is a non-root vertex of the 
graph. If w is not in the graph, then we require that f{v) = f{w). If the root vertex is drawn 
so as to contain the point w, and the root row and column are deleted from L{G), then we 
require that f{w) = 0. We leave it to the reader to check that if / satisfies these boundary 
conditions, then the restriction of / to G is an eigenvector of the matrix with eigenvalue A. 

For certain subgraphs of the square-octagon graph, exact formulas for the number of 
perfect matchings can be found by doing urban renewal to get a weighted version of one 
of the graphs shown below. In each case the eigenvectors of the weighted version can be 
obtained from the eigenvectors of the unweighted version by multiplying by weights 2^~^^ and 
2"*"-' as in § ^. The eigenvalue in the weighted version is obtained from the unweighted 
eigenvalue by adding 1. 

In the rest of this section we give a number of graphs G, their duals G"*", the graphs 
'^(^*)/*)) the eigenvectors of L{G) (possibly with root row and column removed), and 
the corresponding formula for the number of perfect matchings of ?i(f *,/*). When G is 
undirected, we can also enumerate the matchings of 7i(f *, /*) by counting the spanning 
trees of the dual of G. 
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6.1. Temperley's bijection 

Temperley's original bijection involved computing the number of perfect matchings of a 2i—l 
by 2m — 1 subgraph of the square grid with a corner removed. By the main theorem, this 
is the number of spanning trees of an £ x m rectangle. The eigenvectors of the negative 



T — V — 9 — V" 

<) O 6 G- 

<? — — — 0- 



Figure 9: Temperley's original bijection, for the odd-by-odd rectangular region with a corner 
removed. £ = 5, m = 4. 



Laplacian are of the form 



y) = fj,k{x, y) = COS —— cos 



m ' 



where x runs from ^ to £ — | by integer steps and y from ^ to m — ^ by integer steps, j is an 
integer in [0, i— 1] and k is an integer in [0, m — 1]. (The lower left vertex is {x, y) = (|, |) and 
the upper right vertex is {x,y) = (£ — |, m — |.) The eigenvalue of fj^k is 4 — 2 cos — 2 cos 
which is zero when k = £ = 0. The number of spanning trees is 



— 1 I 4-2cos^ - 2cos — 

where the product is taken over all pairs (j, k) with 0<j<£ — l,0<A;<m — 1 except 
(0,0). 

Using panel (b) of Figure ^, we can compute the number of spanning trees in a different 
way. In this case we take the negative Laplacian of the graph in panel (b) and remove a row 
and column corresponding to the outer vertex. The eigenvectors of the resulting graph are 
as follows. 

. vrjx . -Kky 

Jj,k[x, y) = sm — — sm , 

I m 

where this time x runs from to £ by integer steps and y from to m by integer steps, j is an 
integer in [1, ^ — 1] and k an integer in [1, m — 1]. Note that the x, y coordinates correspond 
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to the centers of the faces in the coordinates of panel (a). The corresponding eigenvalue is 
4 — 2 cos ^ — 2 cos — . The number of spanning trees is 



£-1 m-1 

JJ JJ [4 - 2 COS y - 2 COS 

j=l k=l 



m 



In this case there is no zero eigenvalue since we have removed a row and column of the 
negative Laplacian. The equivalence of these two formulas follows from the well-known 
identity for i = U^l\i2 - 2 cos f ). 



6.2. Dimers on an even-by-odd rectangle 



T V T 9 T 

T ? Y T T 

<? — — — — 





Figure 10: The even-by-odd rectangular region, i = 5, m = 4. 



See Figure |T0[ In this case we again remove the row and column from the negative Lapla- 
cian corresponding to the "extended" vertex. The eigenvectors of the negative Laplacian are 
of the form 

7r(2j + l)x irky 

fj,k{^,y) = cos- 



cos ■ 

2i+l m 

where x runs from | to £ + | and y from ^ to m — j is an integer in [0, ^ — 1] and k is an 

-1 m— 1 



integer in [0, m — 1]. The corresponding eigenvalue is 4 — 2 cos ^^^i!^}^ — 2 cos The number 
of spanning trees is 



n-p-r I ^ 7r(2j + 1) ^ ■nk^ 
II 4 - 2 cos ^r^^ — — - 2 cos — 



2£+ 1 



jr=0 /c=0 

Similarly, in panel (b) with the extended vertex removed the eigenvectors are 
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where x runs from to £ and y from to m, j is an integer in [0, £ — 1] and k is an integer 



in [1, m — 1]. The corresponding eigenvalue is 4 — 2 cos ^2^f'/^ ^ ^ 
6.3. Dimers on an even- by-even rectangle 



Trfc 



<> — — — — 



T — 9 — 9 — V — 9 

T T T T T 

-0 — — — 9 — 





Figure 11: The even-by-even rectangular region. £ = 5, m = 4. 



See Figure |Tl]. In panel (a) the eigenvectors of the negative Laplacian after removal of 
the extended vertex are 

7r(2i + l)x 7r{2k + l)y 

fj,k{x, y) = cos — : COS ■ 



2i + l 



2m + 1 



where x runs from | to £ + | and y from ^ to m + |, j is an integer in [0, — 1] and k is an 
integer in [0,m — 1]. The number of spanning trees is 



£-1 m-l 

nn[ 

j=0 k=0 



71(2? + 1) 7r(2k + 1] 

4 - 2 cos ^ . ^ - 2 cos ■ ^ ' 



2i+ 1 



2m + 1 



For eigenvectors in panel (b) replace the cosines with sines, the x range to to i, and 
the y range to m. The formula for the determinant is identical. 



6.4. Dimers on an odd-by-odd rectangle with an extra vertex 



See Figure |T^. In panel (a) the eigenvectors of the negative Laplacian after removal of the 
extended vertex are 

njx . nky 



fj,ki^,y) = COS' 



sm 



m 
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where x runs from ^ to i — ^ and y from to m, j is an integer in [0, ^ — 1] and k is an 
integer in [l,m — 1]. The number of spanning trees is 

TT TT [4 - 2 COS ^ - 2 COS — . 

L £ mi 

j=0 k=l 

The eigenvectors in panel (b) are more comphcated. 
Remark The formula for this region can also be derived by multiplying m (the number 
of ways the extra vertex can be paired) by Temperley's formula for the odd-by-odd region 
with a corner removed, since Temperley's formula still holds with other perimeter vertices 
of the right parity removed instead of the corner. 



6.5. A diamond-shaped region 



See Figure |T^. In panel (a) the eigenvectors of the negative Laplacian after removal of the 
extended vertex are 

7rj{x + y) .^TTk{x — y) 



fj,k{x,y) 



sm 



sm 



2i 2m 

where the origin (x, y) = (0, 0) is in the center bottom part of the extended outer vertex of 
G, X & [—m, i], and ?/ G [0, £ + m], with < x + y < 2i and < y — x < 2m. The indices 
{j, k) G [1, 2£ — 1] X [1, m — 1] U [1, £] X {m}. The number of spanning trees is 

nr. A 
['^-^'''Yi'''2^. 



where {j, k) runs over this range. 
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Figure 13: A diamond-shaped region with £ = 4, m = 3. 
6.6. Another diamond-shaped region 



Figure 14: Another diamond-shaped region, with i = A,m = 3. 

See Figure |1^. In panel (a) the eigenvectors of the negative Laplacian after removal of 
the extended vertex are 

. T^j{x + y) . 7rk{x-y) 
ti k\xy) = sm sm , 

where the origin (x, y) = (0, 0) is at the bottom vertex of G-^, x G [—m, i], and ?/ G [0, £ + m], 
with 0<x + y <2i and < y - x < 2m. The indices (j, k) G [1, 2£ - 1] x [1, m - 1] U [1, £ - 
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1] X {m}. The number of spanning trees is 



4 — 4 cos — cos 

L 2£ 2m 



where (j, fc) runs over this range. Note that this is exactly the formula of the previous section, 
except for the index (j, k) = (£, m) whose eigenvalue is 4. As a consequence the number of 
spanning trees of this diamond is exactly one-fourth the number of spanning trees of the 
previous diamond. This fact was first observed by Stanley (1996) (when i = m), and was 
first proved by Knuth (1997) and Ciucu (1998) (without assuming i = m); a generalization 
was stated and proved by Chow (1997). 

Again we cannot compute the eigenvectors in panel (b). 

Similar formulas hold for two additional diamond-shaped regions with boundaries inter- 
mediate between those shown in Figures O and IH. 



6.7. A triangular region 




o- 
o c>- 

9 — V — Y 
• 9 9 y 

— O A Ch- 

ni) — o — o — o- 




Figure 15: A triangular region, m = 6. 



See Figure |T5|. In panel (a) the eigenvectors of the negative Laplacian after removal of 
the extended vertex are 



7r(2j + l)x 7r(2k + l)y 7r(2fc + l)x 7r(2j + l)y 
sm sm sm sm 



2m + 1 



2m + 1 



2m + 1 



2m + 1 



where the origin [x, y) = (0, 0) is the lower-bottom-most part of the extended external vertex, 
X and y run from to m, with x > y, and j and k are integers in [0, m — 1] with j < k. The 
number of spanning trees is 



n 

0<j<A:<n 



7r(2? + 1) 7r(2fc + 1) 

4 - 2 cos : / - 2 cos ^ ^ 



2m 



2m+l i 
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We don't know how to compute the eigenvectors of the negative Laplacian of the graph 
in panel (b). 

Remark This family of regions has been studied by Mihai Ciucu and Lior Pachter. Ciucu 
(1997) showed combinatorially that the number of domino tilings of the 2n x 2n square 
equals 2" times the square of the number of domino tilings of a particular region so 
one way to prove the claim (not the one given here!) is to hitch a ride on the formula for 
the number of domino tilings of the square — or you could turn this around and derive the 
formula for the number of tilings of the square as a corollary of the above formula. Also, 
Pachter (1997) gave a purely combinatorial proof that the number of domino tilings of Hn 
is odd. A different way to see this uses 2-adic analysis of the factors in the double product 
(Cohn, 1999). 

6.8. Another triangular region, and the quartered Aztec diamond 




See Figure |T^. 

In panel (a) the eigenvectors of the negative Laplacian after removal of the extended 
vertex are 

vrjx Tiky irkx njy 

Jj kix, y) = sin sm sm sm , 

mm mm 

where the origin (x, y) = (0, 0) is the lower-bottom-most part of the extended external vertex, 
X and y run from to m, with x > y, and j and k are integers in [1, m — 1] with j < k. The 
number of spanning trees is 

nr. ^ '^i ^ T^k' 
4-2COS — -2cos — . 
L m m - 

0<j<k<m 
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We don't know how to compute the eigenvectors of the negative Laplacian of the graph 
in panel (b). 

Remark A formula for the number of spanning trees of the "quartered Aztec diamond" 
(G"*" here) had been an open problem for some years before Mihai Ciucu found an expression 
for it (in work that has not yet been written up). Here we effectively get a formula for it by 
doing Fourier analysis on the dual graph. The equivalence of Ciucu's formula and ours does 
not appear to be an effortless identity. 

6.9. A subregion of the hexagonal lattice 

Let uj = e^'^*/'^ and Z[u;] be the lattice of Eisenstein integers. Let G be the directed graph 
whose vertices are each vertex v having an edge directed towards the three vertices 

V + l,v + u,v + u^^, as in Figures |^a and [T^a. 

Here we obtain an exact formula for the number of directed spanning trees on the tri- 
angular region in G whose vertices are 0,m, —muj~^, with "wired" boundary conditions 
(Figure 0). 




Figure 17: A triangular region of the hexagonal lattice. Here n = 5, and m = n — 2 = 3. 

Let n = m + 2. Let A„ be the sublattice of Z[uj] generated by n(l — uj^^) and n{u! — 
uj^'^). A fundamental domain for A„ is given by the hexagon whose vertices are 0,n,n{l — 
uj^-^), —2nuj~^ ,n{uj — uj~^),nuj. This fundamental domain is tiled with six copies of the 
equilateral triangle whose vertices are 0, n, —nuj~^. The graph is embedded in this triangle 
as the set of vertices not touching the boundary. 

Let fa,i3,-y be the function on Z[uj] such that fa,i3,'y{x + yuj + zlj~^) = a^P^'j^; single- 
valuedness implies a/37 = 1- The function /a,/?,^ is an eigenvector of the negative Laplacian 
of G, and has eigenvalue 3 — a — /3 — 7. In order for fa,f3,'r to be a function on the torus 
G/An, the relations (0/7)" = 1 and (P/'y)"' = 1 must also hold. If a is a 3nth root of unity, 
and a/7 is an nth root of unity, then /5/7 = l/(a7^) is also an nth root of unity, so that 
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/a,/3,7 is an eigenvector of the torus G/An with eigenvalue 3 — a — /? — 7. The various fa,i3,yS 
formed in this way are orthogonal on the torus, and since there are Sn^ of them, the same 
as the number of vertices in the torus, they form an eigenbasis of the negative Laplacian of 
the torus. 

Let Li,L2,L^ be the three lines of symmetry of the lattice 7j[w] defined by Li = M, 
L2 = cuM, and L3 = cj^^M. Orthogonal reflection in each of these lines is not only a 
symmetry of the lattice but preserves A„ and the edge directions. Therefore these reflections 
are symmetries of the underlying graph on the torus. The lines Li, L2, L-^ project to the torus, 
cutting it into 6 equilateral triangles. Say that a function on the torus is skew symmetric if 
reflecting it through any of the lines Li, L2, is equivalent to negating the function. Since 
reflection in a line of symmetry Li, L2, or L3 sends the eigenvector /a,/3,7 on the torus to 
the eigenvector f 0,^,(3, f'y,i3,a, or //j,^,^ respectively (in each case two subscripts have been 
transposed), when we express any skew symmetric function as a linear combination of the 
/a,7,/3's, we find that it is a linear combination of the ga,'y,i3^s, which are defined by 

5'o,/3,7 — /a,/3,7 fa,-y,l3 + f'y,ce,l3 f "1,13, a + //3,7,a fp,a,"i- 

ga,i3,'y is an eigenvector of the torus G/An with eigenvalue 3 — a — /3 — 7. If we restrict 
our attention to triples (a,/3, 7) satisfying a < (3 < j under some arbitrary ordering of the 
roots of unity, then the (y'Q-,/3,-y's still span the space of skew symmetric functions; furthermore 
a nontrivial linear relation among them would imply a nontrivial linear relation among the 
/a,/3,7's, SO the 5'a,/3,7's such that a < j3 < 'J form an eigenbasis of the space of skew symmetric 
functions. 

Since the ga,^,"/^^ are zero on the lines of symmetry Li, L2, L^, when we restrict the ga,i3,-y^s, 
to one of the 6 equilateral triangles, they form an eigenbasis of the matrix obtained from the 
negative Laplacian of the triangular subgraph by deleting the vertices on the boundary. Thus 
multiplying their eigenvalues gives the number of spanning trees of the triangular region with 
wired boundary conditions, rooted at the wired boundary: 

n {s-a-p-^y/'. 

a,/3,7 
0/37=1 
o3" = l 

W/3)"=l 
a, P, 7 distinct 

The prime factors of the numbers given by this formula exhibit some rather nice patterns. 
7. Open Problems 

Our work suggests (or might be useful in the solution of) a number of different problems. 

First, there is the natural question of extending the correspondence between spanning 
trees and matchings so that it applies to matchings of more general planar graphs (or perhaps 
even some non-planar ones). An intriguing example is the "12-6-4 lattice" (the infinite graph 
obtained by taking the 1-skeleton of the Archimedean tiling of the plane by dodecagons, 
hexagons, and squares) shown in Figure 0. In a matching, each vertex is paired to another 
vertex via one of three types of edges: a 6/4 edge (bordering a hexagon and square), a 
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Figure 18: A portion of the "12-6-4 lattice" of dodecagons, hexagons, and squares. 



12/4 edge (bordering a dodecagon and square), or a 12/6 edge (bordering a dodecagon 
and hexagon). In random perfect matchings of suitably defined subgraphs (chosen so as to 
minimize the effect of the boundary), the probabilities of these three events are respectively 
1/6 + (19/78)i?, 1/3 + (l/39)i?, and 1/2 - (7/26)i?, where R is given by 



n=0 fc=0 



2n\ frA"^ lA^ 



n \k 132" 



Curiously, these probabilities also turn up in the random spanning trees of a certain directed 
weighted cartesian lattice. In this lattice each leftward edge has weight 25, each upward edge 
has weight 25, each rightward edge has weight 14, and each downward edge has weight 1. The 
probability that in the tree the parent of a vertex is to the left is the same as the probability 
of a 12/4 edge in the 12-6-4 lattice, and the probability that the parent is above the given 
vertex is the same as the probability of a 6/4 edge in the 12-6-4 lattice. These "coincidences" 
suggest a weighted bijective correspondence analogous to the one for the lattice of octagons 
and squares, but we have been unable to find one. 

One side-issue of a number-theoretic nature concerns the formula given at the end of 
§ |6.9| . As was mentioned, there are some interesting patterns governing the prime factors of 
these numbers; for instance, it appears that most of the large prime factors, including the 
very largest, are congruent to 1 mod n. It would be good to know why this is true. 

Turning to entropy, it would also be desirable to have a more general understanding of 
asymptotics. In all the examples considered in § ^ the logarithm of the number of spanning 
trees, divided by the size of the system, tends to a single limit that is independent of the 
shape of the boundary, but is a numerical characteristic of the infinite square grid. One 
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might ask the same question for more general grids. For example, suppose we have some 
locally finite infinite planar graph G that admits an adjacency-preserving action of 1? that 
has only finitely many vertex-orbits and cdgc-orbits. What is the limit of the normalized 
logarithm of the number of spanning trees of large finite subgraphs of G, as the subgraphs 
grow to fill G? It is not too hard to guess the shape of the answer, using the spectrum of 
the negative Laplacian, but justifying this answer would require some care. 

There is also the issue of uniqueness of the measure of maximum entropy. In the case 
of the square grid, Burton and Pemantle (1993) proved that there is a unique translation- 
invariant measure on spanning trees of the infinite square grid that achieves the entropy 
bound. This yields an analogous result on domino tihngs of the plane. It is natural to ask 
the question for lozenge tihngs of the plane; by the theorems of this article, this is closely 
related to the problem of determining whether there is a unique measure of maximum entropy 
for directed spanning trees in the directed triangular lattice. 

The Ising partition function and the number of spanning trees of a graph are both evalua- 
tions of the Tutte polynomial of the graph at special points; see Welsh (1993) for background 
on the Tutte polynomial. Evaluating the Tutte polynomial at most points is #P-complete, 
but at special points the polynomial can be evaluated in polynomial time. For planar graphs 
these special points include the values for the Ising partition function and the number of 
spanning trees. For planar graphs there are bijective connections between Ising systems and 
perfect matchings (Fisher, 1966), and between spanning trees and perfect matchings. Thus 
it is natural to search for connections between perfect matchings and the other special points 
of the Tutte polynomial. 
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